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Developing a theoretical framework for conducting electronic fluids qualitatively distinct from those 
described by Landau's celebrated Fermi liquid theory is of central importance to many outstanding 
problems in condensed matter physics. Perhaps the most important such pursuit is a full microscopic 
characterization of the high-Tc cuprate superconductors, where the so-called "strange metal" behavior 
above Tc near optimal doping is inconsistent with being a traditional Landau Fermi liquid. Indeed, a 
microscopic theory of such a strange metal quantum phase could possibly shed new Hght on the interest- 
ing low-temperature behavior in the pseudogap regime and on the d-wave superconductor itself. Here, 
we present a theory for a specific example of a strange metal, which we term the "d-wave metal." Using 
variational wave functions, gauge theoretic arguments, and ultimately large-scale DMRG calculations, we 
establish compelling evidence that this remarkable quantum phase is the ground state of a reasonable mi- 
croscopic Hamiltonian: the venerable t-J model supplemented with a frustrated electron ring-exchange 
term, which we study extensively here on the two-leg ladder. These findings constitute one of the first 
explicit examples of a genuine non-Fermi Hquid metal existing as the ground state of a reahstic model. 



Over the past several decades, experiments on strongly cor- 
related materials have routinely revealed, in certain parts of 
the phase diagram, conducting liquids with physical prop- 
erties qualitatively inconsistent with Landau's Fermi liq- 
uid theory.^ Examples of these so-called non-Fermi liq- 
uid metal^ include the strange metal phase of the cuprate 
superconductors^ and heavy fermion materials near a quan- 
tum critical point However, such non-Fermi liquid behavior 
has been notoriously challenging to characterize theoretically, 
largely owing to the failure of a weakly interacting quasipar- 
ticle description. It is even ambiguous to define a non-Fermi 
liquid, although possible deviations from Fermi liquid theory 
include, for example, violation of Luttinger's^ famed volume 
theorem, vanishing quasiparticle weight, and/or anomalous 
thermodynamics and transport P'^'CII This theoretical quandary 
is rather unfortunate as it is likely prohibiting a full un- 
derstanding of the mechanism behind high-temperature su- 
perconductivity, as well as stymying theoretically-guided 
searches for new exotic materials. 

Pioneering early theoretical work on the cuprates relied on 
two main premises, '^E^'Cll from which we will be guided but 
not constrained in our pursuit and understanding of a partic- 
ular non-Fermi liquid metal: (1) that the microscopies can be 
described by the square lattice Hubbard model with on-site 
Coulomb repulsion, which at strong coupling reduces in its 
simplest form to the t-J model; and (2) that the physics of 
the system can be faithfully represented by the "slave-boson" 
technique, wherein the physical electron operator is written as 
a product of a slave boson ("chargon"), which carries the elec- 
tronic charge, and a spin- 1/2 fermionic "spinon," which car- 
ries the spin, both strongly coupled to an emergent gauge field. 
However, within the slave-boson formulation, it has been dif- 
ficult to access non-Fermi liquid physics at low temperatures 



because this requires the chargons to be in an uncondensed, 
yet conducting, quantum phase, i.e., some sort of the elusive 
"Bose metal." Early attempts to describe the strange metal in 
this framework treated it as a strictly finite-temperature phe- 
nomenon in which t he sla ve bosons form an uncondensed, but 
classical, Bose fluid,^^a treatment which precludes the pos- 
sibility that the strange metal is a true quantum phase at all. 

In our view, the strange metal should be viewed as a gen- 
uine two-dimensional (2D) quantum phase, which can per- 
haps be unstable to superconducting or pseudogap behav- 
ior. Indeed, recent experimental work on La2-a:Sra;Cu04 has 
shown that when superconductivity is stripped away by high 
magnetic fields, strange metal behavior persists over a wide 
doping range down to extremely low temperatures.^^ Thus, 
the strange metal in the cuprates is quite possibly a true, ex- 
tended, zero-temperature quantum phase.^ 

Inspired by these results and building on our previous work 
which proposed^ and succeeded in realizinj^^'^ a genuine, 
zero-temperature Bose metal, we employ a novel variant of 
the slave-boson approach to construct and analyze an exotic 
2D non-Fermi liquid quantum phase, which we refer to as 
the "d-wave metal" and abbreviate as "d-metal." The d-wave 
metal is modeled by a variational wave function consisting of 
a product of a (i-wave Bose metal wave function^ ^ for the 
chargons and a usual Slater determinant for the spinon. Impor- 
tantly, placing the chargons into the d-wave Bose metal state 
imparts the many-electron wave function with a sign structure 
qualitatively distinct from that of a simple Slater determinant, 
and in particular, imprints strong singlet d-wave two-particle 
correlations. This results in a gapless, conducting quantum 
fluid with an electron momentum distribution function which 
exhibits a critical, singular surface that violates Luttinger's 
volume theorem,^ as well as prominent critical Cooper pairs 



with d-wave character. The d-wave nature of our phase is 
tantaUzingly suggestive of incipient d-wave superconductiv- 
ity and thus of possible relevance to the cuprates. 

Furthermore, tying back into premise (1) above, we pro- 
pose a reasonably simple model Hamiltonian to stabilize the 
(i-metal by augmenting the traditional t-J model with a four- 
site ring-exchange term K. Indeed, such multisite interactions 
are likely necessary to sufficiently frustrate the electron charge 
and qualitatively alter the electronic phase, from a Fermi liq- 
uid to a non-Fermi liquid. Then, owing to the afforded numer- 
ical and analytical tractability provided by the density matrix 
renormalization group (DMRG)^^^^ and bosonization,'^3Eni 
we place the problem on a quasi-one-dimensional (quasi- ID) 
two-leg ladder geometry (see Fig. [T]). In this system, we es- 
tablish several lines of compelling evidence that the d-metal 
phase exists as the quantum ground state of our t-J-K model 
Hamiltonian, and we are able to characterize and understand 
the phase very thoroughly. We believe this study to be one 
of the first unbiased numerical demonstrations of a non-Fermi 
liquid metal as the stable ground state of a local Hamiltonian. 
Finally, in our concluding remarks, we will discuss straight- 
forward extensions of these results to two dimensions, as well 
as comment on their potential relevance to the actual non- 
Fermi liquids observed in experiment. 



GAUGE THEORY AND VARIATIONAL WAVE FUNCTIONS 
FOR THE d-WAVE METAL 

Our theoretical description of the non-Fermi liquid d-metal 
begins by writing the electron operator for site r and spin 
state s = t, i as the product of a bosonic chargon 6(r) and 
fermionic spinon /s(r); that is, Cs(r) = 6(r)/s(r). With 6(r) 
a hard-core boson operator, this construction prohibits doubly 
occupied sites, an assumption we make from here on. The 
physical electron Hilbert space is recovered by implement- 
ing at each site the constraint 6^(r)6(r) = /j(r)/s(r) = 
c|(r)cs(r) = ne(r), which physically means that a given 
site is either empty or contains a chargon and exactly one 
spinon to compose an electron. Theoretically, this is achieved 
by strongly coupling the 6's and /'s via an emergent gauge 
field.E 

Under the natural assumption that the spinons are in a Fermi 
sea state, the behavior of the chargons determines the result- 
ing electronic phase. Condensing the bosonic chargons so that 
(6(r)) 7^ implies Cs(r) ~ /s(r); thus, in this case, the elec- 
tronic phase is that of a Fermi liquid. It then follows that 
in order to describe a non-Fermi liquid conducting quantum 
fluid within this framework, we require that the chargons not 
condense, (6(r)) = 0, yet still conduct. However, access- 
ing such a "Bose metal" phase has proven extremely difficult 
over the years. In recent workj^^^^^'^'we have indeed succeeded 
in realizing a concrete, genuine Bose metal phase, which we 
named the "d-wave Bose liquid" or, equivalently, "d-wave 
Bose metal" (DBM). The DBM is central to our construction 
of the d-wave metal. Specifically, in the DBM, we decom- 
pose the hard-core boson as 6(r) = (ii(r)(i2(r) with the con- 
straint (ij(r)(ii(r) = d2{v)d2{v) = 6'''(r)6(r), where di (^2) 




FIG. 1: Schematic of the t-J-K model Hamiltonian. Top: Picture 
of the full t-J-K model, Eq. ([6|, on the two-leg ladder. In this work, 
we use periodic boundary conditions in the long (x) direction for all 
calculations. Bottom: Action of the ring term Hk, Eq. ([8]), on a 
single plaquette, elucidating its "singlet-rotation" nature. 



are fermionic slave particles ("partons") with anisotropic hop- 
ping patterns', di (d2) is chosen to hop preferentially in the x 
(y) direction. The resulting bosonic phase is a conducting, yet 
uncondensed, quantum fluid, which is precisely the phase into 
which we place the charge sector of the d-metal. That is, for 
the d-metal we take a novel all fermionic decomposition of 
the electron. 



Cs{r) = di{T)d2{r)fs{r), 



(1) 



subject to the constraint 



4(r)di(r) = 4(r)d2(r) = ^/](r)/,(r) = ne(r). (2) 

S 

The resulting theory now includes two gauge fields: one to 
glue together di and d2 to form the chargon and another to 
glue together b and / to form the electron. In the Supplemen- 
tary Information, we give a detailed bosonization analysis of 
this novel gauge theory for the two-leg ladder study on which 
we focus below. 

Guided by the slave-boson construction, one can natu- 
rally construct electronic variational wave functions by taking 
the product of a hard-core bosonic wave function with a 
fermionic wave function f and evaluating them at the same 
coordinates (Gutz wilier projection): 



M{r]h {rt}) = Vg k({Rj) X V^^({rl}, {rt}) 



(3) 



where Vg performs the projection into the physical electronic 
Hilbert space: {RJ = {rj} Ul^^tl- we put the /'s into a 
spin-singlet Fermi sea state with orbitals {kj } (Slater determi- 
nant), i.e., V^/({rJ}, {r^}) = det[e^^^-^I] det[e^^^-^t] = ^js^ 
then we can model both the Fermi liquid metal and the non- 
Fermi liquid (i-metal in a unified way. In both cases, the wave 
functions are straightforward to implement using variational 
Monte Carlo (VMC) methods.^3'^ 

For the Fermi liquid, we put the 6's into a superfluid wave 
function ijj^^ via a typical Jastrow form, so that, schemati- 
cally, i/j^^ = 7^g[V^5^ ^ ^T^' Note that since i/j^^ is a 
positive wave function, the sign structure^^ of ip^^ is iden- 
tical to that of the noninteracting Fermi sea state. In contrast. 
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to model the d-metal, we put the 6's into a Bose metal wave 
function according to the DBM construction of Refs. a2T1 - [24l 

= i^d, ({RJ) X Vrf. ({RJ) = i^F"""^, (4) 

where ipdi (i^d2 ) is a Slater determinant with a Fermi sea com- 
pressed in the x (y) direction. Then, we have 

(5) 

Interestingly, this construction, Eq. ([5]), is actually a time- 
reversal invariant analog of the composite Fermi liquid de- 
scription of the half-filled Landau level,^^ where the d-wave 
Bose metal wave functioiPH plays the role of Laughlin's u = 
1/2 bosonic state. Just as Laughlin's wave function imprints 
a nontrivial complex phase pattern on the Slater determinant, 
the DBM wave function imprints a nontrivial ti-wave sign 
structure, hence our designation "d-wave metal." As we ex- 
plore in detail below, there are many physical signatures asso- 
ciated with putting the chargons into the DBM phase, making 
the d-wave metal dramatically distinguishable from the tradi- 
tional Landau Fermi liquid. 



MICROSCOPIC RING-EXCHANGE MODEL 

The t-J-K model Hamiltonian which we propose to stabi- 
lize the (i-metal phase is given by 

H = Hu^Hk, (6) 
Htj = -t Yl (4c..+H.c.)+J^5^,-5^„ (7) 

i7K = 2K^(5j3524 + H.c.), (8) 
□ 

where (i, j) and □ indicate sums over all nearest-neighbor 
bonds and all elementary plaquettes of the 2D square lattice, 
respectively. In the spirit of the t-J model, we choose to work 
in the subspace of no doubly occupied sites, but for simplicity, 
we do ignore the term — ^n^n^ present in typical definitions 
of the t-J model.El In Eq. ([s]), we have defined a singlet cre- 
ation operator on two sites as Sjj = ^(^J^c^^ — cj^c^^), so 
that Hk can be viewed as a four-site singlet-rotation term (see 
Fig.[T]). For K > 0, the ground state of Hk on a single pla- 
quette with two electrons is a d^^^ -orbital spin-singlet; thus, 
loosely speaking, Hk has a tendency to build in (i-wave cor- 
relations in the system and qualitatively alter the sign structure 
of the electronic ground state. Further arguments for studying 
this model in our search for the (i-metal can be found in the 
Supplementary Information. 

Although seemingly ad hoc for our purposes, we note that 
terms such as our ring term Hk are actually present in careful 
derivations of t- J-type models for Coulomb-interacting elec- 
trons, but are typically ignored largely for theoretical conve- 
nience. In fact, estimating the strength of K, or coefficients 
on related terms, in real materials such as La2-a;Sra;Cu04 is 
an interesting open question. 



TWO-LEG STUDY: DMRG AND VMC 

Unfortunately, as with any interacting fermionic model, 
our t-J-K Hamiltonian suffers from the so-called "fermionic 
sign problem," rendering quantum Monte Carlo calculations 
inapplicable.^^ Also, recently developed tensor network state 
methods^ are as of yet unable to handle critical 2D states 
with a singular surface (e.g., the Fermi liquid and the d-metal) 
due to the anomalously large amount of spatial entanglement 
present.^^ We thus follow the heretofore 

SUCCeSSful2^it24i41j42| 

approach of accessing such phases by studying their quasi- 
ID descendants on ladder geometries, relying heavily on 
large-scale DMRG calculations. In fact, we have already 
established^^ that for two, three, and four legs, the DBM 
phase itself is the stable ground state of a boson ring-exchange 
model analogous to Eq. (|6]). Here, we take the important first 
step of placing the electron ring t-J-K model on the two-leg 
ladder in search of a two-leg descendant of the (i-metal. 

For concreteness, we now consider the model, Eq. ([6]), on 
the two-leg ladder (see Fig. [T]) at a generic electron density 
of p = Ne/{2L^) = 1/3, where TVe = + TV^i is the 
total number of electrons and Lx is the length of our two- 
leg ladder (i.e., the system has x 2 total sites). At this 
density, p = 1/3 < 1/2, on the two-leg ladder, the noninter- 
acting ground state is a spin- singlet wherein electrons of each 
spin partially fill the bonding band (ky = 0), leaving the anti- 
bonding band (ky = tt) empty. Thus, for t K,we expect 



di d2 
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FIG. 2: Picture of the parton bands for the d-metal phase. We 

show orbitals for a 48 x 2 system, showing partially occupied bond- 
ing {ky = 0) and antibonding {ky = tt) bands for di and partially 
occupied bonding bands for d2 and that is, each Slater determi- 
nant in Eq. ^ consists of momentum- space orbitals as depicted here. 
The total electron number is A^e = + — Ndi — Nd2 — 
Nf^ + Nf^ 32, with Nf^ Nf^ 16 so that 5'tot 0; the 
longitudinal boundary conditions are periodic for di and antiperiodic 
for d2 and This is precisely the same d-metal configuration for 
which we display characteristic measurements in Fig. [5] 
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the system to be in a simple one-band metallic state, which 
is a two-leg analog of the Fermi liquid. Formally speaking, 
this phase is a conventional Luttinger liquid with two one- 
dimensional (ID) gapless modes (central charge c = 2). For 
moderate values of ring exchange, K > t, we anticipate the 
unconventional non-Fermi liquid d-metal to be a candidate 
ground state. On the two-leg ladder at this density, the d-metal 
phase has characteristic band filling configurations for the di , 
d2, and /^/^ partons as shown in Fig. [21 di partially fills both 
bonding and antibonding bands, while d2 and only fill 
the bonding band. (The di and ^2 configurations constitute 
the phase denoted "DBL[2,1]" in Ref. 22 ) In a mean-field 
approximation in which the partons do not interact, the sys- 
tem has five ID gapless modes corresponding to the five total 
partially filled bands. However, in the strong-coupling limit 
of the full quasi- ID gauge theory (see the Supplementary In- 
formation for details), two orthonormal linear combinations 
of the original five modes are rendered massive, leaving an 
unconventional Luttinger liquid with c = 3 gapless modes. 

We now provide extensive numerical evidence that this two- 
leg descendant of the d-metal exists as the ground state of 
the t-J-K model over a wide region of the phase diagram. 
We summarize these results in Fig. |3]by presenting the full 
phase diagram in the parameters K/t vs. J/t as obtained by 
DMRG calculations on length Lx = 24 and 48 systems at 
electron density p = 1/3. For small K, we find a conventional 
one-band (spinful) Luttinger liquid phase which is a two-leg 
analog of the "Fermi liquid metal," hence the label in Fig. [3] 
For moderate J and upon increasing K, the system goes into 
the unconventional "non-Fermi liquid d-wave metal" phase, 
which is the main focus of this work. The phase boundaries 
in Fig. [3] all of which represent strong first-order transitions, 
were determined by measuring several standard momentum- 
space correlation functions in the DMRG (see the Supplemen- 
tary Information for details): the electron momentum distri- 
bution function (c^^Cqg), the density-density structure factor 
(^nq(5n_q), and the spin-spin structure factor (Sq • S_q). 

For concreteness, we now focus on the cut along J/t = 2 
in Fig. |3] for a 48 x 2 system with Ne = 32 electrons. We 
take one point deep within the conventional one-band metal 
at K/t = 0.5 and the other point deep within the exotic d- 
metal at K/t = 1.8. First focusing on the former case, we 
show in Fig.|4]DMRG measurements characteristic of the con- 
ventional Luttinger liquid. The ground state is a spin-singlet 
with a sharp singularity in the electron momentum distribution 
function at Qy = and Qx = = ttNc^/Lx = 8 • 27r/48, 
which is a usual Fermi wavevector determined solely from the 
electron density. The density-density and spin-spin structure 
factors at Qy = also exhibit familiar features at Qx =0 
and Qx = 2kF = ■ 27r/48, both characteristic of an or- 
dinary one-band metallic state with gapless charge and spin 
modes. We stress that, even with the constraint of no double- 
occupancy and nonzero K/t = 0.5 and J/t = 2, the interact- 
ing electronic system is still qualitatively very similar to the 
two-leg free Fermi gas; analogously, the 2D Fermi liquid is in 
many ways qualitatively similar to the 2D free Fermi gas. In 
both cases, the main diff erences are basically quantitative and 
are well-understood.'^'^ 




FIG. 3: Phase diagram of the t-J-K electron ring-exchange 
model at electron density p = 1 /3 on the two-leg ladder. In addi- 
tion to the conventional one-band metal ("Fermi liquid metal") and 
exotic "non-Fermi liquid d-wave metal," there are two other realized 
phases. For small J, there is an intermediate phase with fully polar- 
ized electrons (region labeled "FP"). For large K, due to the inher- 
ently attractive nature of ring-exchange interactions !^^^ the system 
generally phase separates along the ladder (region labeled "PS"). 



We turn now to the characteristic point within the d-metal 
phase at J/t = 2 and K/t = 1.8. In Fig.jsj we show a set of 
DMRG measurements at this point, as well as measurements 
corresponding to a variational wave function chosen such that 
its singular features best reproduce the DMRG data (see the 
Supplementary Information for details of our VMC methods). 
The selected d-metal wave function configuration is depicted 
schematically in Fig. |2] Specifically, we have the following 
parton Fermi wavevectors: 2A:^]^ = 21 • 27r/48, 2/c^j^ = 
11 • 27r/48, 2kFd2 = 32 • 27r/48, and 2kFf = 16 • 27r/48. 
The overall agreement between the DMRG and VMC mea- 
surements is very compelling, and we now summarize our un- 
derstanding of these measurements from the perspective of the 
d-metal theory. 

In sharp contrast to the conventional Luttinger liquid, the 
electron momentum distribution function now has singulari- 
ties for both Qy = and qy = tt at a wavevector Qx = = 

i^Fdi ~ ^Fdi]/^- ^^^^ wavevector corresponds to a composite 
electron made from a combination of parton fields consisting 
of a right-moving di parton, a left-moving ^2 parton, and a 
right-moving spinon: d^^^ d2Lf^R' In fact, these "enh anced 
electrons" can be guessed from simple "Amperian rules"'^^^'^ 
in our quasi- ID gauge theory as described in detail in the Sup- 
plementary Information. 

The corresponding density-density and spin-spin structure 
factors, displayed in Fig. |5jb)-(c), also show nontrivial be- 
havior. We expect the density-density structure factor to be 
sensitive to each parton configuration individually and thus 
have singular features at various "2/cf" parton wavevectors 
(see Refs. 21 23 and the Supplementary Information). In the 
DMRG measurements, the most noticeable features are at 
Qy = and Qx = 2/c^j^, which allow us to directly 

read off the realized di parton configuration (see Fig. [2]). The 
lack of these features in the VMC data, as well as the lack of 
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FIG. 4: DMRG measurements in the conventional Luttinger liq- 
uid phase at J/t = 2 and K/t = 0.5. We show (a) the electron mo- 
mentum distribution function, (b) the density-density structure factor, 
and (c) the spin- spin structure factor. The important wave vectors kp 
and 2/cf, as described in the text, are highUghted by vertical dashed- 
dotted lines. 



analogous features at Qx = 2kFd2 in the DMRG data, can be 
understood within our gauge theory framework as presented 
in the Supplementary Information, where we also note that 
our wave function is only a caricature of the full theory. Fi- 
nally, the spin- spin structure factor at Qy =0 not only has 
a familiar, expected feature at Qx = 2k Ff coming from the 
spinon, but also remarkably contains a feature at Qx = 2 Kg 
that can be thought of as a "2/cf" wavevector from the domi- 
nant "electron" in Fig.[5ja). All in all, as we detail thoroughly 
in the Supplementary Information, the DMRG measurements 
are amazingly consistent, even on a fine quantitative level, 
with being in a stable non-Fermi liquid d-metal phase. 

It is important to note that the wavevector Kg depends 
on the interaction strength K/t since the wavevector s k^^-^ 



FIG. 5: DMRG measurements in the unconventional d-metal 
phase at J/t = 2 and K/t = 1.8. We show the same quantities 
as in Fig. |4] Here, we also show the matching VMC measurements 
using a d-metal trial wave function depicted in Fig. [2] 



and k^Fdi ^^^y ^^^^ ^^^S exchange .1^ In Fig. ^ we show at 
J/t = 2 evolution with K/t of the wavevector Kg, i.e., the 
location of the sharp steps in the electron momentum distribu- 
tion function [see Fig.[5ja)], as determined by DMRG. Since 
the momentum- space "volume" enclosed by these singular 
features depends on the interaction K/t and is not simply de- 
termined by the total density of electrons, we may confidently 
say that the (i-metal violates Luttinger's volume theorem.^ In 
fact, the very notion of a single "Fermi surface" is actually am- 
biguous in the d-metal phase. We also show in Fig.|6] for those 
values of K/t at which they are discernible, the wavevec- 

tors 2k p^-^ and 2/c^j^ as identified by features in the DMRG- 
measured density-density structure factor at = [see 
Fig. [5jb)]. For all points, the locations of the identified fea- 



tures satisfy the nontrivial identity Kg = — 2/c^j-^]/4, 
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FIG. 6: Evolution of singular wavevectors in the d-metal phase. 

At fixed J/t — 2 and varying K/t, we show the location of the 
dominant singular wavevector Kg in the electron momentum distri- 
bution function [see Fig. [5|a)], as well as the wavevectors identi- 
fied as 2A;^]-L and 2/^^^-^ in the density-density structure factor [see 
Fig.[5jb)]. These calculations were done with DMRG. 

as predicted by our theory. 

A remarkable property of the d-metal state found in the 
DMRG is that it has prominent critical d-wave Cooper pairs 
residing on the diagonals, as anticipated earlier from the ring 
energetics. We detail these findings in the Supplementary In- 
formation, while here we only mention that such Cooper pair 
correlations have the slowest power law decay among all the 
discussed observables, including the electron Green's func- 
tion. This is in stark contrast with the conventional metal and 
suggests that the d-metal phase has some incipient d-wave su- 
perconductivity in two dimensions. 

As a final piece of "smoking gun" evidence that the real- 
ized DMRG phase is in fact the (i-metal, we have measured 
the number of ID gapless modes, i.e., the effective central 
charge c, via scaling of t he bip artite entanglement entropy ^^ 
in the DMRG and VMCP^^ wave functions. As explained 
above, we expect c = 2 in the conventional Luttinger liq- 
uid and c = 3 in the d-metal. A detailed comparison of the 
DMRG and VMC entropy measurements is presented in the 
Supplementary Information, where the DMRG- VMC agree- 
ment is just as impressive as it is for the more traditional mea- 
surements of Fig. |5] The effective central charge versus Kit 
at J/t = 2 as determined by the DMRG is shown in Fig. [7] 
Indeed, these measurements indicate that c 2 in the con- 
ventional one-band metal, while c 3 in the exotic d-metal. 
Importantly, since 3 > 2, this result demonstrates that our d- 
metal phase is non-perturbative in that it cannot be understood 
as an instability out of the conventional Luttinger liquid. 

DISCUSSION AND OUTLOOK 

In this paper, we have presented exceptionally strong ev- 
idence for stability of a two-leg descendant of our exotic 
strange metal-type phase, the d-wave metal, and we conclude 
here with an outlook on exciting future work. Firstly, it would 
be desirable to march towards two dimensions by studying 
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FIG. 7: Central charge c as a function of interaction K/t, By 
measuring the von Neumann entanglement entropy 5*1 in the DMRG, 
we calculate the effective central charge c at fixed J /t — 2 and vary- 
ing K/t. There is a dramatic jump from c 2 to c 3 at the 
transition, as predicted by our theory. Data for two example points, 
K/t — 0.8 and 1.8, is shown in the inset, where X is the number of 
rungs in each bipartition. (See also the Supplementary Information.) 



systems with more legs. One of the main purposes of Ref.l24l 
was to establish stability of the d-wave Bose metal, the main 
ingredient of the d-metal, on three- and four-leg ladders, and 
this was indeed achieved. Therefore, we do not envision any 
obvious conceptual obstacles in the way of realizing a similar 
result for the d-metal; however, we do anticipate that going 
to more legs will be very challenging numerically due to the 
large amount of spatial entanglement present. In addition to 
more legs, it would also be worthwhile to explore other elec- 
tron densities, in particular those closer to half-filling (p = 1) 
as appropriate to the cuprates. 

With the goal of connecting to experiment, it would also 
be interesting to perform a detailed energetics study of the t- 
J-K model in two dimensions and explore applicability of 
such models to strongly correlated materials. By studying 2D 
variational wave functions based on the d-metal, it should be 
possible to compare physical properties with experimentally 
observed strange metals, such as that in the cuprates. This 
could include various instabilities of the d-metal, e.g., spinon 
pairing as a model of a pseudogap metal or chargon pairing 
as a model of an "orthogonal metal" discussed recently. It 
would be particularly exciting to investigate incipient d-wave 
superconductivity of the cuprate variety — this is rather nat- 
ural given the (i-wave sign structure already inherent in the 
nonsuperconducting parent d-metal. Finally, while we have 
thus far stressed its Luttinger volume violation as a charac- 
teristic non-Fermi liquid property of the (i-metal, we note that 
the 2D phase will also have no Landau quasiparticle as well as 
exhibit non-Fermi-liquid-like thermodynamics and transport. 
Comparing these predictions with properties of real strange 
metals would be an interesting endeavor. In the end, however, 
we would like to stress the conceptual nature of the present 
study, and we hope that our ideas may open up new avenues 
for thinking about non-Fermi liquid electronic fluids. 
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SUPPLEMENTARY INFORMATION 

I. NUMERICAL METHODS AND SUPPORTING DATA 

In this part, we provide some details of the DMRG and 
VMC methods and also show comparison of the Renyi en- 
tropies between the two. We will present additional data anal- 



ysis in Sec. |IIE| after developing a long- wavelength descrip- 
tion of the d-metal phase. 



A. DMRG calculations and observables 

We determine the ground state phase diagram of the t-J-K 
model, Eq. ^ of the main text, by large-scale DMRG calcu- 
lations. We consider square lattice clusters with total number 
of sites Lx X Ly. Here, we study the two-leg ladder system, 
i.e., Ly = 2, and use periodic boundary conditions along the 
X direction. 

Our DMRG calculations generally keep between m = 
5,000 and 20,000 states in each DMRG block. This is found 
to give excellent convergence in the measurements such as 
the ground state energy and various correlation functions de- 
fined below, with small errors which can be neglected safely 
for our sizes up to = 48. The phase boundaries in the 
( J/t, iir/t)-parameter space are determined by extensive scans 
of the derivatives of the ground state energy and by monitor- 
ing the correlation functions. On the other hand, as we will 
describe later, even with such a large m we can converge the 
entanglement entropy only for sizes up to L^^ = 36. 

To characterize the ground state properties of the system, as 
well as properties of the variational wave functions, we calcu- 
late the electron Green's function 



(9) 



where e = Cs with 5 = t, i the electron spin (there is no im- 
plied summation over s). For our electron Green's function 
calculations and analysis, we fix s to one of the two possi- 
ble flavors of spin, say 5 = t, which in the spin- singlet states 



considered in this work, gives the same Green's function as 
the other flavor of spin. The Fourier transform gives us the 
electron momentum distribution function 



(c^gCqg) 



1 ^g.q.(r,-r,)^^t^ 



LxLy 



(10) 



Similarly, we calculate the electron density-density structure 
factor in momentum space 



(JnqJn_q) 



LxLy 



{{ni-p){nj -p)), 



(11) 



where rii = J2s ^Is^^s is the electron number operator at site 
i and p is the electron density. To characterize the magnetic 
properties of the system, we also study the spin structure fac- 
tor 



1 ^giq.(r.-r,)(S,.S.), 



LxLy 



(12) 



where the spin operator is defined as 

1 ^r^ t 

2 2^s,s' ^is^ss'Cis' ' 

Finally, we measure Cooper pair correlations 

Gcooper[r„r,;r;,,r;] = {P[i, j\^P%l]) , 



(13) 



where a Cooper pair operator residing on some nearby sites i 
and j is defined by 



^[hj] = {ci^Cj^ - Ci^Cj^) = P[j, i] . (14) 

We will specifically be interested in diagonal d-wave 
Cooper pairs as detailed in Sec. |IID[ Their correla- 
tions can be accessed by considering the combina- 
tion Gcooper[(^,l),(x + l,2);(x^l),(x' + 1,2)] - 
Gcooper[(^,l),(^ + l,2);(x^2),(x' + 1,1)]; iMs this 



specific quantity which we plot below in Suppl. Fig. 12 



B. Details of the VMC calculations 

The central task of the VMC analysis is to construct a pool 
of variational states for which to collect data that can be com- 
pared with the DMRG results. The determinantal VMC with 
the d-metal wave functions is straightforward: Measurements 
in the VMC simulation of the energy and correlation functions 
described above are all accomplished by averaging over rela- 
tive probability amplitudes corresponding to one or two par- 
ticle hops. These are simply ratios of determinants with one 
or two changed columns, which can be computed efficiently 
using the previously established methods in Ref . 3 1 The one 
exception is the density-density correlator, which merely re- 
quires averaging over the product of density operators. These 
are measured simply by checking for the presence of particles 
on the relevant sites and therefore involve no relative proba- 
bility amplitudes. 



8 



In principle, one can be quite exhaustive and consider all 
possible band fillings of the two bands for each of the four 
partons in the (i-metal wave function, Eq. ([5| of the main 
text, on the 2-leg ladder. Additionally, there exists the free- 
dom to adjust the boundary conditions for each of the par- 
tons so long as the overall conditions are periodic for the elec- 
tron wave function. There are four combinations that accom- 
plish this. Listing the boundary conditions in order of di, d2, 
fs, we have: (1) periodic-periodic-periodic, (2) antiperiodic- 
antiperiodic-periodic, (3) periodic-antiperiodic-antiperiodic, 
(4) antiperiodic-periodic-antiperiodic. Finally, we can intro- 
duce two other variational parameters in the form of exponents 
on the magnitudes of the di and ^2 determinants.^^ That 
is, we can replace V^rfi({Ri}) |V^rfi({Ri})|^i" 
and^,2({Rz}) ^ |V^d2({Rj)|^^-Vd2({Rz}). 

In practice, we do not need to consider all states exhaus- 
tively. From our mean-field understanding of the (i-metal 
phase, we can make two assumptions outright. First, we only 
need to fill the bands with a single contiguous strip in each 
band (centered around = to respect the lattice inversion 
symmetry). Second, we only need to fill the ky = band 
for the spinons. Furthermore, the DMRG results give us three 
important hints. First, in the d-metal region of the phase dia- 
gram, the exact ground states are either spin singlets or have 
a small spin. For the cut J/t = 2, where we focused ex- 
tensively in our analysis, all d-metal states have 5* = or 
5* = 1. Second, the DMRG states show a strong aversion 
to nonzero momentum, which, in the VMC language, corre- 
sponds to symmetric filling of all the bands. Finally, by read- 
ing off the locations of the singular features in the DMRG data 
(specifically, the density-density structure factor), we can de- 
duce the corresponding band fillings for the two bands in the 
di parton Fermi sea. As we increase K/t in the d-metal re- 
gion at fixed J/t, the singular features move such that the di 
band fillings evolve accordingly. For example, for the 48 x 2 
system with Ne = 32 presented in Fig.[6]of the main text, we 



need band fillings ranging from 



(0) 
dl 



(tt) 

^di ~ To accommodate the zero-momentum 
constraint, the up and down spinon numbers oscillate between 



(0) 
dl 



23 and TV, 



(tt) 

dl 



9 to 



N. 



ft 



N- 



fi 



16 for M?^ and TV^^ odd and TV 



(tt) 



'dl 



Nf^ = 15 for Nj^^^ and 7V^^^ even. This corresponds to 
switching between boundary condition type (3) and type (2) 
listed above and total spin S = and 5 = 1, respectively. 

Amazingly, we actually observe this precise shell-filling ef- 
fect in the DMRG data itself. Specifically, in the 2k^^^l = 
N^^^ ' 2ii/Lx data presented in Fig. y of the main text, 
we measure a total spin = when A^^^^ is odd, and 
5* = 1 when N\-,^ is even (since in our calculations — 



dl 



17, 



(0) 



dl 



A/"^^^ is even, A"^^^ and A'^^'^ have the same parity). 
This is very remarkable because the DMRG is producing an 
unbiased ground state of the t-J-K model and, within its in- 
ner workings, has no notion of d and / partons. Thus, a priori 
there is no reason to expect such a precise oscillation of the 
total spin between = and 6* = 1 as we vary the Hamilto- 
nian parameter K/t. Nonetheless, we do indeed observe this 
perspicuous effect, exactly as predicted by our d-metal theory. 



r(0) 



(tt) 



Now going back to the VMC wave functions themselves, 
once we put all of these considerations together the set of vi- 
able band fillings for the four partons is reduced to a very rea- 
sonable number. For example, for the 48 x 2 system, there are 
only eight such "bare" Gutzwiller states across the J/t = 2 
cut. We use the term "bare" when the exponents on the di and 
d2 determinants are set to one. For each of these states, there 
is then the flexibility to adjust these exponents continuously. 
In attempting to match DMRG results, as in Fig.|5]of the main 
text, after choosing the correct bare state, we tune the expo- 
nents until the features visually match as well as possible. Ad- 
mittedly, this process is a bit subjective, especially since gain- 
ing agreement of some singular features comes at the expense 
of poorer agreement of other singular features. However, an 
energetics analysis, with respect to the varied exponents but 
within a fixed bare state, reveals that the chosen states have 
energies that are very close to the rather robust energetic min- 
ima. Performing energy minimization over all possible orbital 
fillings usually gives a state that is off by one or two filled 
orbitals, but again the energies of our chosen states matching 
the DMRG results are not significantly different from the op- 
timal energy. In Fig.|5j the chosen exponents arepi = 0.7 and 
P2 = —0.4. Indeed a negative exponent on the d2 determinant 
is required to achieve a good match with the DMRG. This cor- 
responds to negating the overly strong Pauli repulsion coming 
from all determinantal factors of the wave function. We can 
gain further understanding of the properties of our VMC wave 
functions from the d-metal gauge theory described in Sec. [ll| 
below. 



C. Entanglement entropy results 



The Renyi entanglement entropies are defined by 



1 



Sa{pA) = z ln[Tr(pS)], 

I — a 



(15) 



where pA is the reduced density matrix of a subregion A of the 
lattice. The limit a ^ 1 results in the familiar von Neumann 
entropy 5'i = -Tr {pA Inp^). 

According to conformal field theory, any given Renyi en- 
tropy is expected to scale a^^^^ 



:fCFT 



In 



Lx . ttX 

— sm — — 



(16) 



for a quasi- ID gapless system in the ground state with pe- 
riodic boundary conditions in the x direction. The subsys- 
tem length X represents the number of contiguous rungs con- 
tained in the subsystem, the central charge c is equal to the 
number of gapless modes, and c^^ is a non-universal constant. 
On the two-leg ladder at p = 1/3, we expect c = 2 gap- 
less modes in the conventional Luttinger liquid phase. In con- 
trast, the bosonization approach to the unconventional d-metal 
phase (see Sec.|ll]) predicts c = 3 gapless modes, as there are 
five partially-filled bands but two are rendered massive in the 
strong-coupling limit of the gauge theory. Thus, the central 
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charge as determined by measuring the scaUng of the entan- 
glement entropy provides a crucial diagnostic for the presence 
of the non-Fermi liquid d-metal phase. 

The DMRG has information about the full entanglement 
spectrum and hence can calculate any Sa including the von 
Neumann entropy Si, which is the focus of Fig.[7]of the main 
text. However, calculations of Sa with DMRG are also very 
difficult to converge for highly entangled gapless systems such 
as the d-metal state, and while we were able to well converge 
the von Neumann entropy for Lx = 24 and 36 length systems, 
we were unable to obtain convergence for the Lx = 48 sys- 
tem. Specifically, for the Lx = 48 system at J/t = 2 and 
K/t = 1.8 as focused on in the main text, e.g., in Fig.jsj fit- 
ting to the von Neumann entropy gives c = 2.75 even when 
keeping m = 20,000 states. We estimate that the entropy data 
near the middle of the sample, X Lx/2, is still a few per- 
cent away from full convergence, and, given the sensitivity of 
the fit parameter c on such inaccuracies, we believe a result 
cc^ 3 would be obtained in the limit m ^ oo. 

There are also finite-size uncertainties associated with the 



a scaled-down version of the (i-metal VMC state was used. 



with N^l^ = 15 and N^^^ = 9, exponents pi = 0.7 and 
P2 = -0.4, and Nf^ = TV/; = 12. A linear fit gives c = 2.95 
for the VMC and c = 2.97 for the DMRG, in excellent agree- 
ment with c = 3. It is remarkable that the VMC trial wave 
function is able to reproduce the DMRG data so well and that 
both give the predicted value of c to such high accuracy. 



range of X values used in the fits to Eq. ([16]). For all data 
in Fig. [T] of the main text, we have restricted the fits to data 
-^min = b < X < Lx/2in order to focus on the long-distance 
scaling behavior. Still, at these system sizes, the extracted 
values of c do depend on Xmin, and this by itself produces a 
further uncertainty in c of a few percent. 

In light of the abovementioned uncertainties associated 
with extracting c with the DMRG, and with the goal of further 
bolstering our arguments for the central charge c = 3 within 
the (i-metal on larger systems, we have ventured to compute 
the entanglement entropy in our variational wave functions 
and make concrete comparisons to the DMRG data. Indeed, 
by measuring the expectation value of a "swap operator,"!^ 
VMC can be used to calculatJ^ the Renyi entropies Sa for 
integer a > 2. As such, we now focus on the scaling of ^2 
in both the VMC and DMRG for points within the (i-metal 
phase. 

Let us consider the characteristic (i-metal point examined in 
the main text on the two-leg ladder dXJ/t = 2 and K/t = 1.8, 
with Lx = 48. The scaHng of the VMC and DMRG Renyi en- 
tropies is plotted in Suppl. Fig. [8] There, we observe a strik- 
ing agreement in both the overall magnitudes and detailed fea- 
tures of the two data sets. Indeed, the oscillations in S2, which 
are expected to occur for all but the von Neumann entropy,"^^ 
are in sync with each other. Furthermore, a linear fit gives 
c = 2.95 for the VMC and c = 2.75 for the DMRG, in rea- 
sonably close agreement with the expected value of c = 3 
for the (i-metal phase. For large X, the DMRG entropies lag 
slightly below the VMC entropies: As with the von Neumann 
entropy data discussed above, the DMRG result did not fully 
converge despite having kept m = 20,000 states. (In order to 
attempt to average over the relatively large oscillations in 5*2 
for small X, these quoted values of c used all available data 
Xmin = 2<X <Lx/2m the fits.) 

On smaller lattices, for instance Lx = 36 at J/t = 2 
and K/t = 2 (parameters chosen so that the ground state 
is spin- singlet), the DMRG result is further converged and 
the agreement is even better, as seen in Suppl. Fig. [9] Here, 



II. LONG- WAVELENGTH DESCRIPTION OF THE 
d-METAL PHASE 

In this part, we develop a long-wavelength theory of the d- 
metal phase and provide detailed characterization of the main 
observables, including the electron Green's function, density 
and spin correlations, as well as Cooper pair correlations. We 
illustrate such an analysis with the DMRG and VMC data in 
the (i-metal phase. We conclude with a brief discussion of the 
instabilities of the (i-metal. 



A. Gauge theory description and solution by bosonization 

We start with the decomposition of the electron operator 
into three fermionic partons in Eq. ([T]), as detailed in the main 
text. We assume that the di partons hop preferentially along 
the ladder direction and partially populate both bonding and 
anti-bonding bands, while the d2 partons hop preferentially in 
the transverse direction and populate only the bonding band. 
The densities of the di and d2 partons are the same and equal 
to the total electron density (i.e., including spin). We assume 
that the spinons and populate only the bonding band, 
and the t and | populations are identical so that the state is 
spin-singlet. The parton occupations are illustrated in Fig. [2] 




3.0 



SUPPL. FIG. 8: DMRG and VMC measurements of the scaling of 
the Renyi entropy S2 for Lx = 48 at the characteristic d-metal point 
of J/t = 2 and i^^/t = 1.8. 
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Next, we have the following four-di -chargon term: 



2.5 



SUPPL. FIG. 9: DMRG and VMC measurements of the scaUng of 
the Renyi entropy & for = 36 at J/t = K/t = 2. 



of the main text. The Fermi wavevectors satisfy 



(17) 



where p = Ne/{2Lx) is the total electron density per site. 
When discussing the electron distribution function below, it 
will be convenient to parametrize the di Fermi wavevectors 
as 



knj = 0, TT . 



(18) 



(The same Kg was already used in the main body.) 

We begin with the mean field Hamiltonian in continuum. 
Linearizing near the Fermi points, the kinetic energy density 
reads 

/^kinetic = E E ^^il'^4p^^(-^9x)4p^(19) 

ky^0,7r P=R/L 

+ E P^d2dlp{-id,)d2P (20) 

P=R/L 

+ E E Pvfflp{-id.)fsP ■ (21) 

Here, superscript ky = or tt for di corresponds to bonding 
or anti-bonding respectively; P = R/L = +/— refers to 
right or left moving fields; ^5 is the Fermi velocity for fermion 

"band"6G{4'^; 4"^; ^2; /t; 

Going beyond the kinetic energy, we also consider four- 
fermion interactions, which are required to respect lattice 
symmetries, time reversal, and spin rotation invariance. First, 
we have various density-density interactions in terms of the 
continuum densities p^p for fermion species b. These are 
strictly marginal (harmonic in the bosonized treatment) and 
renormalize band velocities as well as Luttinger parameters. 
Since the final fixed-point theory will have most general form 
consistent with the symmetries, we do not spell out the micro- 
scopic parton density-density terms. 



int,4di 



W 



H.c. 



2w cos 



dl- 



(22) 
(23) 



The last line shows bosonized expression using fields that will 
be defined below; as we will describe later, this non-harmonic 
interaction can destabilize the "relative charge sector" "dl— ". 
We also have the following four-spinon term: 

hintAf = -u (^A^fR^ ■ (/if /l) (24) 

Again, the last line shows bosonized expression using fields 
that will be defined below; this non-harmonic interaction can 
destabilize the "spin sector" "/a". 

The above exhausts residual four-fermion interactions at 



(O/tt) 
Fdl 



On the other 



generic electron densities and generic k 
hand, at commensurate densities we should also include umk- 
lapps. For example, at density p = 1/3 studied in the DMRG 
in the main text, there is a single umklapp term 

^int,umkl = ^umkl ^4i?<^2L ^ fluf^^ + H.C.^ (25) 

= -4^/umki cos(y2i9/^) cos (--^ ^ ^ ~ ^ ' '^^^^ 



If this umklapp term is relevant, it leads to a Mott insulator 
(CDW) of electrons, while it must be irrelevant in the metal- 
lic phase. The absence of any charge ordering in the d-metal 
phase found in the DMRG at p = 1/3 will allow us to crudely 
bound power laws in some correlations and check the overall 
consistency of the d-metal theory, while the umklapp can be 
always rendered inoperative if we step away from the com- 
mensurate density. 

The full theory beyond mean field is a parton-gauge theory. 
Unlike the 2D case, in quasi- ID the parton-gauge theory can 
be solved by bosonization-^^^^^^"^ More specifically, treating 
the gauge field fluctuations on long wavelengths amounts to 
implementing the constraints Eq. ^ for the coarse-grained 
densities along the ladder direction. Since the fixed-point 
Hamiltonian is harmonic in the bosonized fields and the 
coarse-grained densities are linear functions of the bosonized 
fields, implementing the constraints becomes simple. 

Our bosonization treatment is as follows. ES^SI 



bosonize each fermion band 

{4°); 4'^); da; ff, h}) 



species FbR/L (Fb G 



with canonically conjugate boson fields: 

[(pb{x),9b'{x')] = mSbb' 0{x - x') , 



(27) 



(28) 
(29) 
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where Q{x) is the Heaviside step function. Here, we have 
introduced Klein factors, the Majorana fermions {7^5, 7^5/} = 
2^55/ , which assure that the fermion fields with different fla- 
vors anti-commute with one another. The slowly varying 
fermionic densities are simply Pbp = F^pFj^p = dx{P(Pb + 
6>5)/(27r), and hence pb = pbR + pbL = dxOb/i^- 



In the mean field, we start with five modes, 6 



(O/tt) 
dl 



0d2, and 



Of^/l- Including gauge fluctuations is simple in the bosonized 
formulationP^i they effectively implement the constraints 
Eq. ([2]) for the coarse-grained densities and produce the fol- 
lowing pinning of the dual fields (up to fixed constant shifts): 



o(0) 



= % + (30) 

Formally, we perform an orthonormal transformation on the 
fields: 

= 4h (31) 



V2 
1 

71 



(32) 



(33) 



(34) 



(35) 



We perform the same transformation on the ip fields, so the 
new and f fields are again canonically conjugate. Gauge 
field fluctuations render the combinations 9a and 6a in the 
last two lines massive, effectively leading to the pinnings 

and 9n 



ptot — 


2V2 


(%^ 


-% + 


n(0) 






fl(0) , 




ed2) 


Oa = 




+ 




9 A = 










2\/2 





^di 



^di 



in Eq. (30), and only the Of^j 
main 



^di- 



-^ptot 



modes re- 



Note that the definition of 6>ptot is fixed once we define 
the convenient fields 6 fa- and 6di- and require orthogonal- 
ity among these modes and orthogonality to the linear space 
of equations dSOl. On the other hand, the definitions of Oa 



and Oa are somewhat arbitrary (roughly, massiveness of Oa 
corresponds to gluing di and d2 partons to form the bosonic 
chargon b, while massiveness of Oa then glues the b and /); 
we can take any orthonormal linear combination of Oa and Oa 



7ptotJ 



and 



v{djf^-{drOf 



ptot 



(36) 
(37) 



and such a choice does not affect the Of a, 
content of any physical operator below. 

In the fixed-point theory, the "/^" (spin) sector decouples 
from the "dl— " and "ptot" (charge) sectors: 

C = Ci[Ofa]Vfa]gfa\ ^ C[Odi- 
1 

The last line is a generic Lagrangian for a single mode with ve- 
locity V and Luttinger parameter g written in Euclidean space- 
time. In the spin sector, gf^ = 1 is fixed by the condition of 
SU(2) spin invariance. The residual interaction Eq. (24) is 
marginally irrelevant ifu>0. 

On the other hand, in the charge sectors, we have most gen- 
eral quadratic Lagrangian 

C[0^ ^ ' ^ 



'di- 



^ptot 



where we defined 0^ = (Odi-^Optot). In general, A and B 
can be arbitrary positive-definite symmetric matrices. 

To get some feel for the charge sectors, we can consider 
the case where the only interactions between different parton 
species are due to gauge fluctuations enforcing the constraints 
Eq. ( [30| ). For more generality, we allow density-density inter- 
actions within each parton type and encode these in the corre- 
sponding renormalized velocities Vb and Luttinger parameters 
gb. The matrices A and B are readily evaluated to be 



= ^ ^ 



122 



B22 = 



B 
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(38) 



(39) 



(40) 



(41) 



,(42) 



(43) 



As an example, if we take all ^5 = 1 but allow general veloc- 
ities Vb, we can check that the scaling dimension of the inter- 



action Eq. ( 22 ) is greater than 2 and hence this interaction is 
irrelevant. At generic incommensurate electron density, this is 
the only allowed non-linear interaction other than the spin in- 



teraction Eq. ([24]); the latter can be marginally irrelevant, and 
hence the d-metal phase can be stable as a matter of principle. 

For illustrations later in the text, we will consider a 
schematic model where the "dl— " and "ptot" modes also de- 
couple, so 

C = Ci[Ofa;vfa]gfa]^ jC.i[Odi-]Vdi-]gdi-] (44) 
+ jCi [Optot ; Vptot ; ^ptot] . (45) 

This situation arises, e.g., in the above case with gauge-only 



(tt) 



Vdi and gf^ 



interactions, Eqs. (38)-(43), if we further require 

. gdi , where we find 



(0) 
dl 



^dl- 



9ptot 



9dl, 



4 — 



Vdl 

9di 



2^d2 
^d2 



9f 



1 



-1/2 



1 



^dl^dl gd2Vd2 9fVf 



-1/2 



(46) 
(47) 

(48) 



The schematic model holds approximately when the di bond- 
ing and anti-bonding populations are approximately equal. 
Furthermore, we think it is valid for the "bare Gutzwiller" 



wave function, cf. Sec. IB since the projection does not 
know about the intra- species interactions and band velocities; 
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in this case, it is natural to set all bare Luttinger parame- 
ters equal to ^5 = 1 and all velocities equal, and we obtain 
gdi- = Qptot — 1- If we require all bare ^5 = 1 but allow 
general velocities, Eq. (47) would give ^ptot < 1- 



obtained by expanding the sines in the last equation. 



Here we note that the d-metal phase found in the DMRG 
is roughly consistent with the approximation of decoupled 
"(il— " and "ptot" modes with ^ptot significantly larger than 
1 (see our discussion of observables below). Also, as de- 



scribed in Sec. IB our optimal VMC wave functions have 
powers on the determinants and are significantly away from 
the bare Gutzwiller states. We can crudely model the ef- 
fect of adding power on the d2 determinant, V^d2({Ri}) 
|V^d2({Ri})|^'"Vd2({Ri}), by including the Luttinger pa- 
rameter in the bare Lagrangian for the dd2 mode. The 
Slater determinant "11)^2 fills only one band, and a crude guess 
from studies^^ with such ID Jastrow-Luttinger wave func- 
tions is that gd2 = I/P2 for P2 > 0. In particular, 00 
for p2 0. In this limit, ^ptot in Eq. (47 ) can be as large as 2 
even when g^i = gf = 1. Below, we will use values gdi- = 1 
and gptot = 2 to illustrate such a "dressed Gutzwiller" wave 
function. Our VMC measurements with the bare and dressed 
Gutzwiller wave functions support the above conjectures for 
the effective Luttinger parameters in both cases. 

Related to the above discussion, we remark that the DMRG 
does not find any translational symmetry breaking that could 
be driven by the umklapp Eq. (26) at density p = 1/3. In the 
approximation of decoupled "dl— " and "ptot" modes, the ir- 
relevance of the umklapp requires ^ptot > 4/3. This is con- 
sistent with the large effective ^ptot deduced from the DMRG 
correlations and from the optimal VMC wave functions, see 
SecHrElbelow. 



Having discussed the general structure of the d-metal fixed 
point and some approximate models, we now turn to the char- 
acterization of the phase in terms of electron Green's function, 
density and spin correlations, as well as Cooper pair correla- 
tions that are measured in the DMRG. 



B. Electron Green's function 

We begin with the electron Green's function Ge{x^ y), e 
Cs=t/i- I^ the mean field. 



Gf{x,y) 



Gdi{x,y) Gd2{x) Gf{x) , 
sin(/c^J-^^x) sm{kFd2x) sm{kFfx) 



In the first line, we used the fact that the d2 and fs partons pop- 
ulate only the bonding band, hence the ^/-dependence comes 
only from the di parton. In the second line, it is convenient 
to work with bonding or anti-bonding components, ky = 
or TT. The general structure of the oscillating contributions is 



Gq (X, ky 



2A 



(ky) ' 



a^^=± sign{x)\x\ 



(49) 



(50) 



ae'^y Ke + (a + 2/3 + 7)7rp , (51) 



(k ) 

where we allowed general amplitudes A^,^^ and scaling di- 

(k ) (k ) 

mensions A^^^. Contribution at the wavevector Q^^^ comes 
from the following combination of long-wavelength parton 
fields. 



where a,/3,7 = R/L +/-; 5 =t / |; and 



e 



(fey) 

Q!/37,s 



^iky 

72 

Q,g«fej, 



V2 

a- P 



Odl- 



ptot ^ ^ 

a + 2/3 + 7^ 



75 



2V2 

-a - 2/3 + 37 
2V6 



ptot -r ^ 



(52) 

(53) 
(54) 
(55) 
Ofa (56) 
(57) 



The last line contains Oa and Oa that are pinned upon includ- 
ing gauge fluctuations, and from here on we often drop these 
fields when they only modify overall complex phases of the 
operator expressions. For ease of reference. Table |I] lists dif- 
ferent cases where we fix a = +. We will consider these 
cases and entries in the table after some more general discus- 
sion. We are primarily interested in the scaling dimensions 

A (^y) 

First, we note that the spin sector "/c^" is decoupled from 
the rest and has gfcr = 1; the corresponding fields con- 
tribute 1/4 to the scaling dimension of the electron corre- 

(k ) 

lation for any wavevector Q^,^^. Thus, the scaling dimen- 
sion depends only on the "dl— " and "ptot" content of both 

(k ^ (k ) 

Kp-y^s and 6^7^ , fields. We expect that the main difference 
among the cases comes from the 6>ptot content, since its coef- 
ficient (a + 2/3 + 7)/(2V2) can vary in magnitude from for 
a = 7 = — /3 to a/2 for a = 7 = /3; this content is listed in 
one of the columns in Table |T| Note, however, that in general 
the relative signs of the "dl— " and "ptot" components are 
also important and affect the scaling dimensions, thus making 
the scaling dimensions also depend on ky ; this is indeed what 
we find from fitting the DMRG data (see Sec. |II E| ). It is only 
in the case when the "dl— " and "ptot" decouple that the scal- 
ing dimensions do not depend on the signs of the coefficients. 

Second, by using the general result A[e* > 
^ I XI j ^j^j I valid for any canonically conjugate set (pj , 6j and 
any coefficients aj , bj , we can obtain an exact general bound 



Aky) 



> 



1 |2a + 2/3 + 7| 



(58) 
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B 


I 


O^^)^^ Ea J51I 


t^ptot coeii. 

C.+2/3+7 




a/37 ' ' ^^^^ moQei 

Hprniinlpfl — " anH "/ot.nt," 


«^ 5 Dale vjulz wilier 

— 1 r/^+^-t- — 1 

ydi— — -L? i/ptot — -L 


(fc ) 

1 dressed Gutzwiller 

n^-i — 1 n^i-^i- — ^ 
ydi— — -L? i/ptot — ^ 


+ 




+ 







> 1/2 


A+-+ 


1 


3/4 


+ 








1 

V2 


> 1/2 


+^ptot/8 


9/8 


1 


+ 


+ 






1 


> 1 


+^ptot/8 


9/8 


1 


+ 


+ 


+ 


47V p + e'^y Ke 


72 


> 3/2 


+^ptot/2 


3/2 


7/4 



TABLE I: Analysis of different contributions, Eq. p2l, to the electron operator. We list the wavevector Q^^^^ ; the coefficient of Optot in 

(k ) (k ) (k ) 

, which strongly affects the scaling dimension A^^^^ ; exact lower bound on A^^^^ ; the scaling dimension A^^^^ in the approximation 

ot decoupled "dl— " and "ptot" modes, with A+ ^ given in Eq. ( [60| ; specialization to gdi- = Qptot = 1 appropriate for the bare Gutzwiller 

wave function; and specialization to gdi- = 1, Qptot = 2 appropriate for the dressed Gutzwiller wave function that roughly captures the 
DMRG structure factors (see text for details). Note that in general the scaling dimensions also depend on ky. 



which is also listed in Table HI 

Third, we can calculate the scaling dimensions in the ap- 
proximation that the "dl— " and "ptot" modes decouple: 



fluctuations pin the fields 6a and 6a and hence can reduce the 



A+-+ 



h H 7^ 9ptot 



9di- 



32 



(59) 
(60) 



This illustrates the preceding discussion of the dependence 
of the scaling dimensions on the 6>ptot content. By setting 
9di- = Qptot = 1, we obtain scaling dimensions in the bare 
Gutzwiller wave function also listed in the table. On the other 
hand, by setting gdi- = 1 and ^ptot = 2, we obtain scaling 
dimensions in the dressed Gutzwiller wave function (tenta- 
tively corresponding to powers on the determinants pi = 1 
and p2 = 0); these dimensions are listed in the last column 
in Table [11 Power law fits for the electron Green's function 
measured in the VMC are consistent with these predictions. 

Let us now consider different entries in Table Jl We can 
roughly understand the trends in the scaling dimensions by 
appealing to so-calle d ''Am perean" rules for interactions me- 
diated by gauge fields ■^'^^ Parallel currents attract and hence 
processes containing such currents are enhanced, while anti- 
parallel currents repel and hence such processes are sup- 
pressed. The splitting of the electron into three partons leads 
to two gauge fields. We again use the picture where we first 
break the electron into the spinon / and chargon 6, and then 
break the chargon b into the partons di and d2 . We can then 
think of the two gauge fields as follows. The first gauge field 
works to glue the di and ^2 together to form the chargon; the 
two partons carry opposite gauge charges with respect to this 
gauge field, hence processes that contain the di and ^2 mov- 
ing in opposite (same) directions are enhanced (suppressed). 
The enhanced combinations produce the main features in the 
boson distribution function in the so-called DBL[2,1] phase 

of Ref. 22 at wavevectors {kFd2 — k^pdh ^y) i^arked in the 
top panel in Suppl. Fig.[TO| The second gauge field works to 
glue the chargons and spinon to form the electron; the di/2 
and / carry opposite gauge charges with respect to this gauge 
field, and hence processes that contain di/2 and / moving in 
opposite (same) directions are enhanced (suppressed). 

The precise mathematics in (1 + 1)1) is that the gauge field 



fluctuating content in Eq. ( [55] ), particularly for the Amperean- 
enhanced combinations. The first two rows in Table U cor- 
respond to oppositely moving di and d2 partons. From the 
exact lower bound, we see that these entries can be potentially 
more enhanced than the other two rows. In the first row, the 
spinon moves parallel to the di and anti-parallel to the ^2, 
while in the second row the situation is interchanged; both sit- 
uations produce the same exact lower bound, but the first one 
is apparently more enhanced in the approximate model with 
decoupled "dl— " and "ptot" modes. The third row in Table[l| 
has the di and ^2 partons moving in the same direction, but 
the corresponding Amperean suppression is somewhat com- 
pensated by the spinon moving in the opposite direction to 
the two; in the approximate model with decoupled "dl— " and 
"ptot", the second and third row have the same scaling di- 
mension. The last row in Table [l| has all partons moving in the 
same direction; this combination is suppressed by all gauge 
field fluctuations, and we see that it has the largest scaling di- 
mension, which is always larger than the mean field value of 
3/2. 

We conclude with a simple understanding of the main land- 
scapes in the electron distribution function, as measured in the 
DMRG and VMC in Fig. |5];a) of the main text. This discus- 
sion also re-iterates the origin of the wavevectors in the first 
two rows in Table|l] We can take the chargon distribution func- 
tion from the earlier DBL[2,1] study. This has dominant fea- 
tures at wavevectors obtained by combining oppositely mov- 
ing di and d2 partons. For simplicity, we approximate the 
dominant features by ^-functions as shown in the top panel 
in Suppl. Fig. [To] For the spinons /, we take a simple step 



distribution shown in the middle panel in Suppl. Fig. [TO] The 
electron distribution function is a convolution of the chargon 
and spinon distribution functions. The result is shown in the 
bottom panel in Suppl. Fig. [TO] The inner steps arise from 



combinations where / is moving parallel to the di and anti- 
parallel to the d2 and correspond precisely to the first row in 
Table |l| while the outer steps arise from combinations where 
/ is moving anti-parallel to the di and parallel to the ^2 and 
correspond to the second row in Table |l] We can see from Ta- 
ble |l| that the inner singularities are likely stronger than the 
outer ones, particularly if ^ptot is large, and this is consistent 
with the DMRG findings in the main text. 
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SUPPL. FIG. 10: The electronic d-metal is obtained by putting 
the chargons into the DBL[2,1] (Bose-metal) state of Ref. 22 and 
spinons into the Fermi sea (spin Uquid) state. Top panel: Crude ap- 
proximation to the chargon distribution function where we replace 
the peaks in Fig. 8(a) of Ref. 22 with delta-functions; the wavevec- 

(k ) h 

tors are kFd2 — kp^i — Tvp — ^Ke and arise from the oppo- 
sitely moving di and d2 partons. Middle panel: Spinon distribution 
function. Bottom panel: Electron distribution function obtained as 
a convolution of the chargon and spinon distribution functions. The 
wavevectors of the step singularities correspond to the first two rows 
in Table |l| Note that this sketch does not contain potentially impor- 
tant wavevectors listed in the third row in TableH 



have the right symmetry properties: 



^/2 



ir][ sm{V2(pdi-)e ^ 



_u(-) 
Fdl '^Fd 



— "iL "iL 



(61) 
(62) 
(63) 

(64) 
(65) 

(66) 
(67) 



In the second line, ± refers to e*^^. Note also that at 
p = 1/3 we have 2kFd2 = —^kpf, and the product 
^(2/cFd2,o)^(2/ci./,o) gives the umklapp term Eq. ([26]); however, 
at any other density these are distinct wavevectors. Combina- 



tions similar to those in Eqs. ( 63 ) and ( 65 ) but with minus sign 
between the terms contribute to a current operator rather than 
the density fluctuation and are not spelled out here. 

We can similarly obtain fermionic bilinear contributions to 
the electron spin operator, e.g.. 



lJsR 



- sm 



ptot 

e V2 . (68) 



In the bosonized expressions above, we have omitted 6a and 
Oa as these are pinned by the gauge field fluctuations; the dis- 
cussion of the power law correlations in the d-metal phase 
only depends on the displayed fluctuating field content. 

Finally, we have bilinears carrying zero momentum con- 
tributing to both charge and spin, e.g.. 



s,P 



(69) 



such contributions have scaling dimension 1. 

Let us first note some general relations among scaling di- 
mensions for various density and spin correlations: 



C. Density and spin correlations 



A[n. (0) 

y'^Fdl 



1 

2 " 
> 1 . 



, (70) 
(71) 



We now discuss the electron density and spin correlations. 
Erom the various ways of writing the electron number in terms 
of partons, cf. Eq. ([2]), we can immediately obtain fermionic 
bilinear contributions to the electron number operator that 



Thus, the scaling dimension of (o) _,(7r) ^ is always 

y^Fdl ^Fdl^^) 

larger than the mean field value of 1. 

Next, we can get some quantitative feel by using the ap- 
proximation of decoupled "(il— " and "ptot" modes; the seal- 
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ing dimensions in this approximation are: 



A\n 



Qptot 



ddl- Qptot 



Qdi- < 1. This is consistent with the condition that the inter- 
action /iint,4di in Eq. (22) is irrelevant. 



A[n. (0) (.) J = -- 



1 , Qpto 



(72) 
(73) 
(74) 



1 / 1 



^[Vi^-41,^)] = 2^^+^^^-J ' ^^^^ 
A[n(2,,,,o)] = A[5(2,,,,o)] = ^ + ^ • (76) 

As before, we obtain results for the bare Gutzwiller wave 
function by setting gdi- = Qptot = 1- In this case, the 
smallest scaling dimension is A[n(2/c^^2,o)] = \ correspond- 
ing to slow \x\~^ power law. We also have A[n (fcy) ] = 



Furthermore, inspired by the DMRG observation of a sin- 
gularity in the spin structure factor at wavevectors (±2Ke, 0), 
we have also considered four-fermion contributions to the spin 
operator, since, in the strongly coupled theory, these can be 
comparable in prominence to the fermionic bilinears. We have 
identified two interesting terms, which can be constructed by 
combining the already exhibited bilinears. 



5, 



ASf. 



• (0) (tt) 

ir]l ^r]l ' cos 



,^)] = A[n(2/e^^,o)] = A[5(2fe^^,o)] = i cor- 
power law, which is also enhanced over 



: sin(V2l9/^) sm{V2ipdi-) . 



(78) 

(79) 
(80) 



-5/4 



responding to 

the mean field The enhancement of the correlations Note that the first term can be also constructed by combin- 

comes from the reduction of the fluctuating content, w hich is 
the (1+1)D realization of the "Amperean" enhancemen t-^ | 44 | 45 J 
of correlations for processes that contain parallel gauge charge 



.(^)t 



mg c"_^' ' and 



iky) 



and is a kind of "electron 2k f'' from the 



currents. On the other hand, A[n. 



,(0) 



J = 1 and is 



dominant feature observed in the electron distribution func- 
tion. Similarly, the second term can be constructed by com- 



bining, e.g., c^^^^ and cJ^"^ Writing out such electron bilinears 
gives six-parton terms, which turn out to be equivalent to the 
above four-parton terms upon considering the pinning of Oa 
and Oa- In the last equation, we have also carefully kept track 
of these pinned fields, since the very presence or absence of 
this contribution can depend on the pinning values. Both the 
DMRG and VMC have a visible feature in the spin structure 
factor at the wavevector (0, tt), cf. Fig.jsj^c) of the main text, 
and hence we conjecture that the pinning of the fields 6a and 
Oa is such as to give nonzero cosine for their combination ex- 
hibited above. This knowledge will be useful when discussing 
Cooper pair correlations in the next section (note that the pre- 
cise pinning values of the 6a and Oa fields were not important 
in the observables discussed earlier). 

In the same spirit, we can also construct four-parton contri- 
^jj^ butions to the electron density, n(2Ke,o) and n(o,7r) • We expect 
the former to be present generically with properties similar to 



not enhanced over the mean field, since the particle and hole 
partons are on the same side of the Fermi sea and create anti- 
parallel currents. We performed measurements in sample bare 
Gutzwiller wave functions and verified the dominant 
and power laws for different oscillating components 

in the density structure factor. 

We remark that the d-metal realized in the DMRG appears 
to have larger ^ptot > 1 and hence ^(2/cFd2,o) can be sup- 
pressed compared to n (ky) In fact, from our earlier 



discussion of the umklapp term Eq. ( 26 ) at density 1/3 and 
absence of charge order in the d-metal phase found in the 
DMRG, we can conclude that A[/iumki] > 2, and hence have 
an exact bound: 



4 2 2 

A[n(2fc^,,,o)] = g A[/iumki] -9^3- 



In the approximation of decoupled "dl— " and "ptot" modes, 
this corresponds to ^ptot > 4/3, and in this case ^(2feFd2,o) 
already has larger scaling dimension than n (ky) (if we 

,0) 

assume g^i- = 1). Let us also quote the numbers for our 
earlier dressed Gutzwiller example where we set gdi- = 1 
and gptot = 2: We get A[n(2fe^d2,o)] = 1 that is larger 

t^^^^[^(24'^i^,\0)] = = ^[^(2^./,0)] = 

^[S{2kFf,o)] = 3/4, SO the latter singularities are more 
strong. The DMRG indeed finds visible features in the density 

structure factors at wavevectors {2k^pdi ' which we track to 
identify the di bonding/antibonding orbital populations. On 
the other hand, at density 1/3, 2/cf/ = —2kFd2 = 27r/3 and 
these wavevectors are not visible in the density structure fac- 
tors in much of the data. In the above approximation, we 
can make the (2/c^J-^\o) to be dominant if we also assume 



(2Ke,0)- 



On the other hand, the expression for the latter con- 



tains sine of the Oa and Oa combination in Eq. (79), and we 
conjecture that the pinning is such that the n(o,7r) term van- 
ishes. This is consistent with the absence of any feature in the 
DMRG and VMC density structure factor at (0, tt), compare 
Figs, [gb) and (c). 

From the bosonized expressions, we can see that the scaling 
dimension of is related to that of the allowed interaction 



^int,4di,Eq. (22): 



A[h 



int,4diJ 



(81) 



Since the stability of the d-metal requires A[hint,Mi] > 2, we 
conclude that A[S^q ^^^] > 1. This is indeed what we found 
from the DMRG measurements. In the approximation with 
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decoupled "dl— " and "ptot" modes, we find 

2^ 2 ' 
1 1 



^[^f2Ke,0)] 



(82) 
(83) 



To summarize, many such quafitative and semi-quantitative 
considerations of the observables, relations among exponents, 
and stabiHty to perturbations are internally consistent, giving 
us more confidence that the phase found in the DMRG is in- 
deed the electronic d-metal. 



D. Cooper pair correlation 

To complete the discussion of observables, we also consider 
electron Cooper pair operators. Besides being prominent ob- 



servables themselves, these are needed, e.g., in Sec. II F for 
the detailed description of the phases proximate to the (i-metal 
when the electron operator gets gapped. We consider only 
spin singlets and can start with microscopic Cooper operators 



defined as in Eq. (14) on some nearby pair of sites. We can 
write each electron operator in terms of the partons and obtain 
contributions containing six continuum parton fields (two for 
each species di,d2, and /). We examined all such six-fermion 
terms, and below we present most interesting ones bearing in 
mind the application to proximate phases. Rather than doing 
direct expansion, we can use symmetry arguments to identify 
which microscopic Cooper pairs receive particular contribu- 
tions. 

The most important such six-fermion terms are 

VlnL = d2Rd2Ldfld[]lf^LhL + {R^ L) (84) 



(0,7r) 



(0) (tt) 



\/3 



V, 



(0,7r) 



d2Rd2Ldfld^^^f^LfiL -{R^L) 



(85) 



(0) (tt) 

r]\ 'rj\ 'riorums 



2{V20a - eg) 



i2\/20p 



These terms carry momentum (0, tt), i.e., they are translation- 
ally invariant along the ladder and are odd under the leg in- 
terchange. They are defined to be invariant under time re- 
versal, as is appropriate for the singlet Cooper pairs. The 
"(e)" ["(o)"] combination is even [odd] under mirror (x, y) 
{—x^y)\ the numerical constant is sine [cosine] of the par- 
ticular combination of the pinned fields Oa and Oa- We ex- 
pect that one or the other combination is nonzero, but which 
one depends on the detailed pinning specifying the d-metal 
phase (see below). From the identified symmetry properties, 
we can see that the combination 7^(o ^r) contributes to leg- 
bond Cooper pairs that are anti- symmetric under the leg in- 
terchange: 



P[(x, 1), {x + 1, 1)] - P[(x, 2), {x + 1, 2)] ^ {x) - 



On the other hand, the combination contributes to diag- 

onal Cooper pairs anti- symmetric under the leg interchange: 

P[(x, 1), (x + 1, 2)] - P[(x, 2), (x + 1, 1)] ^ V\^^^ (x) + . . . . 

(87) 

The 7^(0, tt) terms have the smallest content of fluctuating fields 
in the d-metal — only the (j)ptot part that is necessary to en- 
code the electrical charge of the Cooper pairs. In the approx- 
imation of decoupled "(il— " and "ptot" modes, the scaling 
dimension is 2/ ^ptot (equal to 2 in the bare Gutzwiller and 
1 in the dressed Gutzwiller wave functions). By measuring 
whether the leg-bond or diagonal Cooper pair shows power 
law behavior, we can constrain the appropriate pinning of the 
fields 6a and 9a- As we will present in Sec. HE both the 
DMRG and VMC find that it is the diagonal Cooper pairs that 
have dominant power law correlations in real space. Hence 
we conclude that the pinning is such that V\'^^^ is nonzero 



ie) 

while V) ' X is zero, which is also consistent with the pres- 



(0,7r) 



ence of 5'(o,7r), Eq. (80), and the absence of n(o,7r) features 



in the DMRG and VMC. The specific diagonal Cooper pairs 



Eq. (87) are natural in the model with electronic ring ex- 
change: e.g., they arise when solving the ring Hamiltonian 
for two electrons on a single placket. They can be viewed as 
having a d-wave character, which is one of the motivations for 
naming our non-Fermi-liquid phase as "d-wave metal." 
We also mention the following six-fermion terms 



(0,0) 



d2Rd2Ldf^^d^^^^ {f^nf^L - hnftL) (88) 



; ( 2 \/2 </)pt ot ± V2 0d 1 _ ) 



cos(V2%,) , (89) 



where ± refers to e'^^y . Such terms carry zero momentum, 
i.e., they are translationally invariant along the ladder and are 
even under the leg interchange. They are also even under 
the mirror (x^y) {—x^y) and contribute, e.g., to a rung 
Cooper pair P[(x, 1), (x, 2)] (as well as leg-bond or diagonal 
Cooper pairs symmetric under the leg interchange). In some 
loose sense, they can be viewed as "s-wave" Cooper pairs, 
while the 7^(o,7r) ones are "d-wave", the precise distinction ly- 
ing in the transformation properties under the discrete ladder 
symmetries. Because of the additional fluctuating field con- 
tent, we expect the 5-wave ones to have larger scaling dimen- 
sion [equal to 2/^ptot + l/(2^(ii-) + 1/2 in the decoupled 
"dl— " and "ptot" approximation] and to be less visible in the 
d-metal than the d-wave ones. However, the s-wave can be- 
come comparably prominent in a spin gap phase discussed in 
Sec.lIlFl 



(86) 



E. Sample of the DMRG power law fits 

Guided by the long-wavelength description of the (i-metal, 
we performed detailed fits for the power laws in various cor- 
relations at dominant wavevectors. Here we highlight some 
results at the DMRG point J/t = 2, K/t = 1.8, presented 
in Fig.|5]in the main text. To remind the readers, the best d- 
metal candidate for this 48 x 2 ladder with 32 electrons has 
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2k^^^ - 



11 • 27r/48, cf. Fig. [2] of the 



21 • 27r/48, 2k^^j^ 

main text. 

The main features in the electron distribution function oc- 
cur at the wavevectors Hsted in the first three rows of Table U 
Ke = 57r/48, 27rp + Kg = 377r/48, and 27rp - Kg = 277r/48, 
for either ky = or ir. We fit the electron Green's function to 
an expression 



Ge{x,ky) = ^ 



a [(L^/tt) sin(7rx/L^)]^^" 



(90) 



The dominant oscillation is at the wavevector Kg, which 
is readily recognized in the electron momentum distribution 
function in Fig. [5|a) of the main text. As illustrated in 
Suppl. Fig. [TT] including only this wavevector already cap- 
tures the overall real-space dependence and gives 2A^1^ ^ 

1.3 and 2A^1^ ^1.5 for the bonding and antibonding elec- 
trons respectively. Fits including the other two wavevectors 
(not shown) capture also finer features without affecting much 
the estimates of A^^^\ We thus confirm the general ex- 
pectation that the exponents can be different in the bonding 
and antibonding electron distribution functions. Supplemen- 



tary Figure 11 also shows the electron Green's function mea- 
sured in the VMC wave function. Similar fits in this case (not 



shown) give roughly equal values 2 A 



(0) 

+-+ 



2A 



(tt) 

+-+ 



1.25, 



which is also expected in the approximate gauge theory treat- 
ment with decoupled "dl— " and "ptot". Our trial wave func- 
tions and the approximate treatment of the gauge theory are 
not qualitatively accurate in this respect. The structure of the 
gauge theory is qualitatively accurate — it is only that we do 
not know numerical values of the couplings in the full theory 
with coupled "dl— " and "ptot" modes. Nevertheless, we see 
that the approximate treatment provides a reasonable semi- 
analytical guide. 

Turning to the spin and density correlations, we fit these 
to an expression similar to Eq. ( [9Q| , but with cosines in- 
stead of sines. It is simple to fit the spin correlations at 
ky = TT, since there is only one feature, Eq. (80). We es- 
timate 2A[S^Q ^^^] ^ 2.2, which satisfies the d-metal stabil- 
ity requirement discussed after Eq. dSTl) and gives us an es- 
timate gdi- ^ 0.85 in the approximation with decoupled 
"dl— " and "ptot". On the other hand, the spin correlations 
at ky = have three important wavevectors: 0, 27r/3, and 
2 Kg = 107r/48. The power law for the zero-momentum com- 
ponent is fixed at Fitting power law decays at the other 
two wavevectors allows us to estimate ^ptot ^ 3.5 — 4 and 
gdi- consistent with the previous estimate. The large value 
of 9ptot explains the weakness of the feature at 27r/3. As dis- 
cussed in Sec. |IB[ our matching VMC state has small nega- 
tive power p2 = —0.4 on the ^2 determinant, which can in- 
deed reproduce such large ^ptot • On the other hand, the value 
of Qdi- < 1 implies that the component at 2 Kg has power 
law decay slower than in real space and singularity in the 
structure factor that is stronger than slope change; the singu- 
larity at this wavevector can be noticed already in the DMRG 
data in Fig. [5jc) of the main text, while the VMC has harder 
time reproducing this. 
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SUPPL. FIG. 11: Real-space DMRG and VMC electron Green's 
functions G'e(x) in the = 48 ladder at J/t = 2, K/t = 1.8, cor- 
responding to the momentum distribution functions shown in Fig. |5] 
of the main text. Thin solid lines show power law fits of the DMRG 
data to Eq. ( [90| restricted to a single oscilating component with 
wavevector Kg. We show data for all distances 2 < x < Lx/2 
to bring out the oscillation more clearly, while the fit itself is done 
over 4 < X < Lx/2 to pick up the long-distance behavior. 



It is also simple to fit the density correlations at ky = 



TT, where we have two wavevectors k 



and k 



(0) 



2A[n. (0) u) 



(0) 
Fdl 



k 



M 

Fdl 



Fdl 



107r/48. The power law fits 



2.3 and 2A[n., (o) 



27r/3 

give 
2.1. 



The latter is consistent with the general bound in Eq. (71 
and with g^i^ value close to 1, while the former is consis 
tent with ^ptot ~ 4. On the other hand, the density correla 



tions at 



have many wavevectors: 0, 2/c^^^\ 



27r/3, and 



2 Kg. Based on the previous estimates, we expect that all these 
components have similar scaling dimension of order 1. We 
can indeed get nice fits, but because of the many parameters 
it is difficult to give accurate individual exponents; neverthe- 
less, all are consistent with the previous estimates of gdi- and 

^ptot • 

Finally, we also measured the diagonal d-wave Cooper pair 
correlations in the DMRG and VMC. For the 48 x 2 sample 
above, we found it difficult to fully converge the pair correla- 
tions in the DMRG (the electron Green's function shown ear- 
lier converged more readily). Nevertheless, we observed very 
clearly a slow power law decay in this Cooper channel, fitting 
roughly - in the DMRG and - in the VMC. 

For illustration, in Suppl. Fig. [12] we show the measurements 
in the smaller 36 x 2 sample at J/t = 2 and K/t = 2, where 
the DMRG data is better converged, with the remaining uncer- 
tainty less than 5% (the same sample was used in the Renyi 
entropy comparisons in Suppl. Fig. [9]). We see good match be- 
tween the DMRG and VMC results, with power law fits giving 
roughly - in the DMRG and - in the VMC. 



Comparing with the predictions in Sec. |IID[ the exponents are 
again consistent with ^ptot value close to 4. All other Cooper 
pair correlators, including also the leg-bond d-wave Eq. ([86]), 
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SUPPL. FIG. 12: DMRG and VMC diagonal d-wave Cooper pair 
correlations [see expression below Eq. |T4] )] in the system of length 
= 36 at J/t = 2 and K/t = 2 (same as in Suppl. Fig.|9]). We 
show the DMRG data for x < 15, where the convergence error is 
estimated to be less than 5% having kept m = 8,500 states in the 
DMRG. Thin line shows fit of the DMRG data to the functional form 
^ [{Lx/tv) sm{7vx/Lx)]~^, withp ^ 1.1. In the inset, we illustrate 
the computed correlation function by showing two diagonal d-wave 
(dxy) Cooper pairs, one being created and the other being destroyed, 
separated by a distance x. 



decay much faster. This observation constrains the pinning 
values of the Oa and Oa fields in our fixed-point theory of the 
d-metal, as discussed with regards to which of the two Cooper 
are present. Furthermore, note that 



pairs Eq. ( [85] ) or Eq. 
the measured power law decay is comparable or even slower 
than that in the electron Green's function, so the diagonal d- 
wave Cooper pair operator is one of the most prominent ob- 
servables in the d-metal phase on the two-leg ladder. It is this 
observation that leads us to speculate about incipient d-wave 
superconductivity in this strange metal phase. 

To conclude, many detailed properties of the various mea- 
sured correlations match well with those expected in the d- 
metal theory. This provides strong support for the identifica- 
tion of the phase found in the DMRG as the d-metal. 



1) If the interaction /iint,4/ is relevant, it pins Of a- and gaps 
the spin. Assuming /iint,4di remains irrelevant, we have two 
gapless modes "(il— " and "ptot" and the full theory is similar 
to the DBL[2,1] theory in Ref. 22 A prominent density ob- 



servable is n 



(2fcF/,0) 



whose wavevector (27rp, 0) 
is determined by the electron density. Loosely speaking, this 
phase can be thought of as a two-leg analog of a "pseudogap," 
i.e., a nonsuperconducting quantum fluid with a spin gap. 

2) If the interaction hint,4di relevant, it pins (fdi-- As- 
suming hint, 4 f remains irrelevant, we have two gapless modes 
"/cr" and "ptot", with one Luttinger parameter in the latter 
sector. In this case, the spin correlations remain gapless. De- 
pending on the precise pinning of (fdi- determined by the 
sign of in Eq. ([22]), we get dominant density correlations 
~ sm{^/2(pdi-)e^^p^°^^^ or current correla- 

cos{V^(pdi-)e'^^p^°^^^', the wavevec- 



tions J. (0) 



tor is (27rp, tt) and is distinct from the case 1). This phase 
can be viewed as a phase where the bosonic chargons form 
a paired-boson state,^^ while the spin remains gapless, and 
is akin to a non-Fermi liquid termed "orthogonal metal" dis- 
cussed recently in Ref. (50l 

3) Finally, if both /iint,4di and /iint,4/ are relevant, we have 
pinning of the fields Ofcr and (fdi- as in the cases 1) and 2) 
above. Only one gapless mode "ptot" remains. We have pres- 
ence of both (27rp, 0) and (27rp, tt) wavevectors in the density 
(or appropriate current depending on the sign of w) correla- 
tion with the same power law. Furthermore, the Cooper pair 
7^(0,0)' Eq. (88 ), at wavevector (0, 0) now becomes as promi- 
nent as the A^,7r) at (0, tt). This phase is the two-leg analog 
of a superconductor. 

Since the DMRG did not observe strong density correla- 
tions at the wavevector 27rp along the ladder or any other sig- 
natures of instability, we conclude that our d-metal phase is 
stable in the t-J-K model. In future work, it would be inter- 
esting to modify the model to further explore the above prox- 
imate phases. 



III. MOTIVATION FOR THE t-J-K MODEL TO REALIZE 
ELECTRONIC d-METAL 



F. Instabilities of the d-metal 

With all observables at hand, we can now consider phases 
that would be obtained out of the d-metal when one or both 



of the interactions hint,Adi and /iint,4/, Eqs. (22) and (24) re- 
spectively, become relevant. Since the electron operator gets 
gapped out, a detailed description requires consideration of 
the Cooper pair operators. In the d-metal phase and in all 
phases below, the Cooper pair 7^(o,7r) ^ 5*^^"^^*°* discussed 
in Sec. |IID| is always prominent, and only one more Cooper 
pair will come to prominence in one of the phases. On the 
other hand, clear distinctions between possible phases are al- 
ready apparent with the density and spin observables. There 
are three cases: 



The purpose of this section is to give some analytical intu- 
ition why the electronic ring Hamiltonian likes the particular 
d-metal phase. Of course, we need direct numerical studies 
to determine what actually happens in the particular Hamilto- 
nian, as done in the main text. 

We can gain some feeling for the energetics by doing a slave 
particle mean field calculation. For convenience, we write the 
electronic ring term as the following electron pair-hopping, 

^1234 = {c\^C2o^ (4^^4/3) + (cLc4a) (4^^2/3^ • (91) 



Here and below, summation over repeated spin indices is im- 
plied. Assuming decoupled di, d2, and / partons with trans- 
lationally invariant hopping mean field, we can calculate the 
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expectation value of the ring term on an elementary square as 

{Rr,r+x,r+x+y,r+y)mi = 2 {\Xf,x\^ + \Xf,y\^) ^ (92) 
X {\Xdl,x\' - \Xdl,y\') {\Xd2,x\' - \Xd2,y\') , (93) 

where Xda,fi = {dl{r)da{r + /i))mf, and x/,^ is evaluated 
similarly for one spinon species (we assume identical hop- 
ping for the t and | spinons). The plus sign between \xx\'^ 
and P for the /-spinons comes from adding the two terms 
in Eq. ( [9T] ) and summing over the spin labels, while we have 
the minus sign for the d-partons from the Fermi statistics. If 
di (d2) partons hop preferentially in the x (respectively y) 
direction, then |xdi,xP > \Xdi,y\'^ (respectively \xd2,x\'^ < 
\Xd2,y\'^)- In this case, the expectation value of the ring term is 
negative, so for the positive K in our model we obtain low ring 
energy; increasing the anisotropy in the di and d2 hoppings 
lowers the ring energy, while the optimal anisotropy is de- 
termined from the competition with the electronic kinetic and 
spin exchange energies. This is our crude mean field argument 
why the ring terms prefer the particular fractionalized state 
and why the di and d2 partons develop strong anisotropics for 
large K. 

We can also provide a more constructive motivation for the 
electronic ring model as a candidate for realizing the d-metal 



phase. This "reverse engineering" argument starts with an ef- 
fective lattice gauge theory for the d-metal phase with partons 
di, (^2, and /. To simplify further, let us go back to splitting 
the electron = bfa into bosonic chargon b and fermionic 
spinon and seek a system where the chargons are in a "d- 
wave Bose metal' (DBM) phase introduced in Ref . |2T1 while 
the spinons are in a spinon Fermi sea state. The latter can be 
obtained by postulating simple hopping for the spinons. On 
the other hand, to put the chargons into the DBM state, we 
follow prior works proposing such a phase in frustrated boson 
models with ring exchanges and realizing it in two-leg, three- 
leg, and four-leg ladder studies. "^^^ Note that this "mean field 
Hamiltonian" for the chargons is actually strongly coupled; it 
was motivated in Ref. [2T] as a candidate that would produce 
the DBM phase where we further split b into di and ^2 partons 
with anisotropic hopping, precisely as we want in the (i-metal 
of electrons. Here by starting with the ring Hamiltonian for 
the 6's we avoid duplicating efforts in Refs. l2TH24l 

Going beyond the "mean field," we need to include a com- 
pact U(l) gauge field that attempts to glue the b and fa partons 
to produce the physical electron.^ Denoting this gauge field 
on the (oriented) lattice links {rr') as Arr' and the conjugate 
integer- valued electric field as Err' , the lattice gauge theory 
Hamiltonian for the d-metal phase has the structure 



H = El, - K^cos(V xA)-Y^ {tiiU-'^^'-' fUr'a + H.c.) - ^ (ttt'e'^-^'&J&r' + H.c.) 

{rr') □ {rr') {rr') 

+ i^^^) ^ J^g^A.,.+.g^A.+,+,,.+, ^ g^A.,.+, g^A.+,+,,.+, ^ blbr+^bl^.^^br+y + H.C. 



(94) 



The Hamiltonian is supplemented by the Gauss law constraint 

(V • ^) (r) = blbr - fUra • (95) 

In the gauge theory Hamiltonian, t^^J and t^^^, are hopping 
amplitudes for the spinon and chargon respectively. The b 
and / partons are coupled to the gauge field A with opposite 
gauge charges, so that the physical electron operator is gauge 
neutral. We have also generalized the boson ring term with 



coupling K^^"^ to a gauge-invariant form that also respects the 
square lattice symmetries. The h and k, terms are standard for 
the lattice gauge field dynamics. 

For large h, the electric fields are pinned at Err' = and 
the partons are confined. We can identify the sector Err' = 
with the physical Hilbert space of electrons. Starting with this 
limit and working perturbatively in t^^\ t^^\ and K^^\ we 
obtain an effective Hamiltonian for electrons that contains the 
following terms: 



E 



^ra^'^'c^ H.C. 



j2 



{rr') 



h 



Sr ' Sr' 



□ 



r-\-x ,r-\-x-\-y ,r-\-y 



H.c. 



. . . .(96) 



The electron hopping is obtained from second-order processes spinons. The electron ring term is obtained from third-order 
hopping both chargon and spinon, while the spin-spin interac- processes involving hopping spinons on opposite edges of a 
tion is obtained from second-order processes exchanging only square and chargon ring exchange on the square [for simplic- 
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ity, we have used isotropic spinon hopping t;^;^^ = t;^;^^ = 

t^^^]. Some other terms are also generated at the same orders 
in but are only indicated with dots; these include terms 
that correlate the electron densities and similarly modify the 
already present terms, but are not qualitatively new. 

The terms exhibited above give precisely the electronic ring 
model studied in the main text, where we allow ourselves to 
vary independently the electron hopping and ring-exchange 
amplitudes, as well as the antiferromagnetic exchange cou- 
pling. In the end, it is the detailed numerical study that estab- 
lishes the phase diagram of this Hamiltonian. 

We conclude by mentioning that the d-wave correlations 
present in our two-leg (i-wave metal are of the dxy variety. 



These are built in by taking the ring term to operate on elemen- 
tary plaquettes of the square lattice: (r, r + x,r + f + ^,r + ^). 
Looking forward, both when going to more legs with the 
DMRG and when studying two dimensions directly with the 
VMC, it will be interesting to consider a c?-wave metal of 
the dx'2-y2 variety, a phase which is potentially more rele- 
vant to the cuprates. A promising model to realize this phase 
would include, instead of that considered in this work, a ring- 
exchange term operating on all plaquettes (r, r + f — ^, r + 
2x,r -\- X -\- y). Addressing the applicability of such mod- 
els to real cuprate materials, as well investigating incipient 
(ia^2_^2-wave superconductivity out of the putative dx2-y2- 
wave metal, are very exciting future problems. 
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